Universal Constraints on Conformal Operator Dimensions 



Vyacheslav S. Rychkov" and Alessandro Vichi^ 

" Scuola Normale Superiore and INFN, Pisa, Italy 
^ Institut de Theorie des Phenomenes Physiques, EPFL, Lausanne, Switzerland 



Abstract 

We continue the study of model-independent constraints on the unitary Conformal Field 
Theories in 4-Dimensions, initiated in arXiv:0807.0004, Our main result is an improved 
upper bound on the dimension A of the leading scalar operator appearing in the OPE of two 
identical scalars of dimension d: 

X = 1 + Oa + . . . 

In the interval 1 < d < 1.7 this universal bound takes the form 

A < 2 + 0.7(d - 1)^/2 + 2.l{d - 1) + 0.43((i - 1)^/^. 

The proof is based on prime principles of CFT: unitarity, crossing symmetry, OPE, and con- 
formal block decomposition. We also discuss possible applications to particle phenomenology 
and, via a 2-D analogue, to string theory. 



1 Introduction and formulation of the problem 



Our knowledge about non-supersymmetric Conformal Field Theories (CFTs) in four dimensions 
(4D) is still quite incomplete. Suffices it to say that not a single nontrivial example is known which 
would be solvable to the same extent as, say, the 2D Ising model. However, we do not doubt that 
CFTs must be ubiquitous. For example, non-supersymmetric gauge theories with A^^^ colors and 
Nf flavors are widely believed to have "conformal windows" in which the theory has a conformal 
fixed point in the IR, with evidence from large Nc analysis [1] , supersymmetric analogues [2] , and 
lattice simulations [3]. Since these fixed points are typically strongly coupled, we do not have 
much control over them. In this situation particularly important are general, model-independent 
properties. 

One example of such a property is the famous unitarity bound [1] on the dimension A of a 
spin / conformal primary operator Oa,z 

A>1 (/ = 0), (1.1) 
A>/ + 2 (/>1). 

These bounds are derived by imposing that the two point function (00) have a positive spectral 
density. 

As is well known, 3-point functions in CFT are fixed by conformal symmetry up to a few arbi- 
trary constants (Operator Product Expansion (OPE) coefficients). The next nontrivial constraint 
thus appears at the 4-point function level, and is known as the conformal bootstrap equation. It 
says that OPE apphed in direct and crossed channel should give the same result (see Fig. [T]). 

The bootstrap equation goes back to the early days of CFT [5]. However, until recently, not 
much useful general information has been extracted from il|^ All spins and dimensions can apriori 
enter the bootstrap on equal footing, and this seems to lead to unsurmountable difficulties. 




Figure 1: The conformal bootstrap equation. The thick red line denotes a conformal 
block, summing up exchanges of a primary operator O and all its descendants. 



Recently, however, tangible progress in the analysis of bootstrap equations was achieved in 
[7]. Namely, it was found that, in unitary theories, the functions entering the bootstrap equations 



^Here we quote only the case of symmetric traceless tensor operators. 

^Except in 2D, in theories with finitely many primary fields and in the Liouville theory We will comment 
on the 2D case in Sections 4.1 and |5] below. 
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Figure 2: The bound f{d) is the smallest Amm for which a positivity property exists. 



(conformal blocks) satisfy certain positivity properties which lead to general necessary conditions 
for the existence of solutions. 

The concrete problem considered in [7], and which we will continue to discuss here, was as 
follows. In an arbitrary unitary CFT a Hermitean scalar primary 0^ of dimension d was singled 
out. The conformal bootstrap equation for its 4-point function {(j)d4>d(pd4>d) was studied under the 
sole assumption that all scalars in the OPE (pa x 4>d have dimension above a certain number, call 

it Ajnin ■ 

0d X 0d = 1 + (Scalars of dimension > Amin) + (Higher spins) . (1-2) 
It was shown that the conformal bootstrap does not allow for a solution unless 

A^in < fid) , (1.3) 

where f{d) is a certain continuous function, computed numerically. We stress that this conclusion 
was reached without making any assumptions about dimensions or spins of other operators ap- 
pearing in the OPE, beyond those implied by the unitarity bounds. Nor any assumptions about 
the OPE coefficients were made (apart from their reality, which is again implied by unitarity). 

In other words, in any unitary 4D CFT, the OPE of any scalar primary (pd must contain at 
least one scalar field Oa with dimension not larger than f{d). 

Incidentally, the function f{d) was found to satisfy /(I) = 2, which is quite natural since d = 1 
corresponds to the free field whose OPE contains the operator :0^: of dimension 2. 



What makes the result like (1.3) possible? The basic reason is that, in any theory, crossing 



symmetry relation of Fig. [T] cannot be satisfied term by term, but only by cancellations among 



various terms. The guaranteed presence of the unit operator in the OPE (1.2) creates a certain 
"crossing symmetry deficit", which has to be balanced by other fields. The idea is to show that 
this cannot happen unless at least one scalar of sufficiently low dimension is present. 

Technically, the method of [7] consists of 3 steps (see Section [3] for a detailed review): 

1. We Taylor-expand the conformal bootstrap equation near the "self-dual point" configuration 
having equal conformal cross-ratios u = v. The expansion is truncated to a certain finite 
order A^. 

2. We systematically search for positivity properties satisfied by linear combinations of Taylor 



coefficients of the conformal blocks, for fields appearing in the RHS of the OPE (1.2). A 
found positivity property implies that the "crossing symmetry deficit" cannot be balanced 
and rules out a CFT with a given d and Amm- 
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3. For fixed d, the bound f{d) is then computed as the point separating those Amin for which 
a positivity property exists, from those ones for which it does not (Fig. |2|. 

The nature of the method is such that increasing can make the bound only stronger. The 
optimal bound should in principle be recoverable in the limit ^ oo. In practice the value of 
is determined by the available computer resources and algorithmic efficiency. The best bound 
found in [7], plotted in Fig. [3| corresponds to = 6. 

m, f=f6 
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Figure 3: The bound fe{d) ^2 + 1.79\/d— l + 2.9{d— 1), I <d < 1.35, corresponding 
to N = 6, reproduced from 




The purpose of this paper is to present an improvement of the bound (1.3) obtained by using 
the method of [7j with larger values of A^, up to A^ = 18. The new results are interesting in two 
ways. First, pure numerical improvement turns out to be significant. Second, A^ = 18 happens 
to be large enough so that we start observing saturation of the bound. So we believe our current 
results are close to the optimal ones achievable with this method. 

The paper is organized as follows. In Section |2] we review the conformal bootstrap equations. 
In Section |3] we review the connection of the bound ( |1.3[ ).with positivity properties satisfied by the 
conformal block expansion coefficients. In Section |4] we present and discuss our results. We also 
mention accompanying results which we obtain for an analogous problem in 2D. In Section |5] we 
propose several future applications and extensions of our method, with emphasis on connections 
to phenomenology and string theory. In Section [6] we summarize and conclude. In Appendix |A] 
we collect some details about our numerical algorithms. In Appendix [B] we include the tables on 
which plots in Section |4] are based. 



2 Review of conformal bootstrap 

We will review the conformal bootstrap equation in its simplest form — as applied to the 4-point 
function of identical scalars (0000). We largely follow [7j, where a more detailed discussion and 
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Figure 4: The operator product expansion of 0(x)0(O) converges for this configuration. 



references can be found. 



2.1 Conformal block decomposition 

Let = 0rf be a Hermitean scalar primarjj^ operator. The operator product expansion (OPE) 
(j) X (f) contains, in general, infinitely many primary fields of arbitrary high spins and dimensions]^ 



4>{x)m 

Here 



\x 



2d 



E 



=2n 



\X\ 



'Ki{x)-O^M + 



Ki{x) 



\x\ 



(2.1) 



I = 2nhy Bose symmetry; 

A > 1 (A > / + 2) for / = (/ > 2) by the unitarity bounds (p^; 



• The ■■• stands for contributions of descendants of the primary Oa,« (i-e. its derivatives). 
These contributions are fixed by conformal symmetry; 

• The OPE coefficients ca,; are real (see Appendix A of [7j). 

We assume that the OPE converges in the following weak sense: it gives a convergent power 
series expansion for any (2 + n)-point function 

{ct>{x)mMyi)---Myn)) 

provided that |x| < i.e. 0(x) is closer to the origin than any other local field insertion (see 
Fig. |4]). This assumption can be justified by using radial quantization ([8], Sect. 2.9), and checked 
explicitly in free field theory. For rigorous mathematical results about OPE convergence see [H]. 

■^The field is called primary if it transforms homogeneously under the 41? conformal group. 
^If there are several primaries with the same A,Z, they have to be all included in this sum with independent 
coefficients. 
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Figure 5: The auxiliary z coordinate. The conformal blocks are regular outside the cut 
denoted by the zigzag line. 



The OPE (2.1) can be used to obtain conformal block decomposition of the 4-point function 



{<jy{x,) <jy{x2) <f>{x,) <f>{x,)) = ' (2-2) 



0:34 



g[u,v) 



1 + X^PA,/^A,i(w,^) , PA,i = cii>0, (2.3) 



where u = ^^2^34/ (^13^24)5 = ^14^23/ (^13^24) the conformal cross-ratios. This representation 
is obtained by using the OPE in the 12 and 34 channels. The conformal blocks gA,i{u,v) sum up 
the contributions of the primary Oa,i and all its descendants. Their explicit expressions were 
found by Dolan and Osborn [10]: 



zz 



9A,i{u,v) = — [kA+i{z)kA-i-2{z) -{z^ z)] , 

Zi z z 

kp{x) = x^'\F, (/5/2, /?/2, /3; x) , (2.4) 
u = zz, V = {1 — z){l — z) . 



Notice the judicious introduction of the auxiliary variables z and z. When the theory is formu- 
lated in the Euclidean space, these variables are complex-conjugates of each other. To understand 
their meaning, it is convenient to use the conformal group freedom to send 0:4 ^ 00 and to put 
the other three points in a plane, as in Fig. [5j Then it's easy to show that 

z=^+X + iY, z = z* , (2.5) 

where (X, F) are the coordinates of x^ in the plane, chosen so that X = F = corresponds to x^ 
halfway between x\ and X3. This "self-dual" configuration, for which u = v, will play an important 
role below. We can see that the z variable is a natural extension of the usual complex coordinate 
of the 2D GET to the 4D case. 

According to the above discussion, the OPE is expected to converge for \z\ < 1. Conformal 
block decomposition is a partial resummation of the OPE and thus also converges at least in this 
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range. In fact, below we will only use convergence around the self-dual point z = 1/2. However, 
conformal blocks, as given by (2.4), are regular (real-analytic) in a larger region, namely in the 



2;-plane with the (1, +oo) cut along the real axis (see Fig. |5]). The conformal block decomposition 
is thus expected to converge in this larger region. One can check that this indeed happens in the 
free scalar theory. 

One can intuitively understand the reason for this extended region of regularity. The condition 
for the OPE convergence, as stated above, does not treat the points x and symmetrically. On 
the other hand, the conformal blocks are completely symmetric in xi ^ X2 and so must be 
the condition for their regularity. The appropriate condition is as follows: the conformal block 
decomposition in the 12-34 channel is regular and convergent if there is a sphere separating the 
points Xi^2 from the points x^^^. For the configuration of Fig. [5| such a sphere exists as long as X2 
is away from the cut. 



2.2 Conformal bootstrap and the sum rule 



The 4-point function in (2.2) must be symmetric under the interchange of any two Xi, and its 



conformal block decomposition (2.3) has to respect this symmetry. The symmetry with respect to 
Xi X2 or X3 X4^ is already built in, since only even spins are exchanged [lUj . On the contrary, 
the symmetry with respect to xi ^ X3 gives a condition 



v'^g{u,v) = u'^giv,u) 



(2.6) 



which is not automatically satisfied for g{u,v) given by (2.3). This nontrivial constraint on di- 
mensions, spins, and OPE coefficients of all operators appearing in the OPE x is known as the 
conformal bootstrap equation. Physically it means that OPE applied in 12-34 and 14-23 channels 
should give the same result (Fig. [T|. 

In the 2;-plane of Section 2.1, the LHS of (2.6) has a cut along (1, -l-oo), while the RHS has a 
cut along (—00,0). Thus, if (2.6) is satisfied, the cuts have to cancel, and the resulting g{u,v) is 
real analytic everywhere except for 2; = 0, 1. 

In [7], we found it useful to rewrite ( |2.6[ ) by separating the unit operator contribution, which 
gives 

u''-v'' = yPA,i [Aa,/(w,^) - WgAAv.u)] . (2.7) 



The LHS of this equation is the "crossing symmetry deficit" created by the presence of the unit 
operator in the OPE. This deficit has to be balanced by contributions of the other fields in the 
RHS. 

In practice it is convenient to normalize (2.7) by dividing both sides by u'^ — v'^. The resulting 
sum rule takes the form: 



1 



d,A,l 



iX,Y) 



Y.Pa,iF,,aAX.Y) , 

v'^gA,i{u,v) - u'^gA,i{v,u) 



{21 
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The "F-functions" -F^ a,; are real and regular in the full 2;-plane cut along (— oo,0) U (l,+oo). In 
particular, the 0/0 behavior at the self-dual point 2; = 1/2 is regular. 

All F-functions vanish near the points z = and z = 1. Thus the sum rule can never be satisfied 
near these points if only finitely many terms are present in the RHS. The OPEs containing finitely 
many primaries are ruled out. 



3 Positivity argument 

The main idea of [7j was very simple, and can be described as follows. Suppose that for a given 



spectrum of operator dimensions and spins {A,/} the sum rule (2.8), viewed as an equation for 
the coefficients pA,i > 0, has no solution. Then of course such a spectrum would be ruled out. 

Any concrete realization of this idea needs a practical criterium to show that there is no 
solution. For a prototypical example of such a criterium, imagine that a certain derivative, e.g. 



dx (see (2.5)), when applied to every F^^a,/ and evaluated at a certain point, is strictly positive 



("positivity property"). Since the same derivative applied to the LHS of (2.8) gives identically 
zero, a solution where all coefficients pA,i are non-negative would clearly be impossible. We refer 
to this simple reasoning as the "positivity argument" . 

One can imagine more general criteria using different differential operators, and applying them 
at different points. In [7j, we found it convenient to apply differential operators precisely at the 
self-dual point z = l/2, X = Y = Q. One can show that the F-functions are even with respect to 
this point both in the X and Y directions: 

F(x,r) = F(x,-r) = F(-x,y). 

Thus, all odd-order derivatives vanish, and a sufficiently general differential operator ("linear 
functional") takes the form: 

A[F]= ^n.,nd^dlF\x=Y=0. (3.1) 

m,n even 
2<m+n<Af 

where is some fixed finite number, and Xm,n are fixed real coefficients]^ Notice the exclusion of 
the constant term m = n = 0, in order to have A[l] = 0. 

Assume that for certain fixed d and Amin, we manage to find a linear functional of this form 
such that ( "positivity property") 

A[Frf,A,«] > for all A > A^i„ (/ = 0) (3.2) 
and for all A > / + 2 (/ = 2, 4, 6 . . .) . 



^In [7j, we analytically continued to the Minkowski space by Wick-rotating Y —> iT. In this picture z and z 
are both real and independent, and conformal blocks are real regular functions in the region < z, z < 1. For our 
purposes Minkowski and Euclidean pictures are exactly equivalent. In particular, derivatives of F-functions in Y 
and T are trivially proportional to each other. 
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Moreover, assume that all but a finite number of these inequalities are actually strict: A.[F] > 0. 
Then the sum rule cannot be satisfied, and such a spectrum, corresponding to a putative OPE 



(1.2), is ruled out. 

The proof uses the above "positivity argument". Since A[l] = 0, the positivity property 
imphes that only those primaries for which A[F] = would be allowed to appear in the RHS of 
the sum rule with nonzero coefficients. By assumption, there are at most a finite number of such 



primaries. However, as noted in Section |2.2[ finitely many terms can never satisfy the sum rule 
globally, because of the behavior near z = 0, 1. Q.E.D. 

While the above formal reasoning is quite sufficient to understand our results, in [7] the sum 
rule was also given an alternative interpretation in terms of convex geometry. In this more vi- 
sual picture, linear combinations of F-functions with arbitrary positive coefficients form a convex 
cone in the space of two-variable functions. One can consider the full function space or its finite- 



dimensional subspace corresponding to Taylor-expanding up to order A^. Positivity property (3.2) 
means that there is a hyperplane separating the function 1 from the convex cone. Thus it im- 
plies that the sum rule cannot be satisfied. The converse is "almost true", modulo questions of 
convergence. 

Clearly, the language of linear functionals provides an equivalent, dual formulation of the 
problem. This formulation is also especially convenient from the point of view of checking our 
results independently. It's not so important how we find the functionals. As long as we publish 



the functional coefficients \m,n-, anyone can verify that the inequalities (3.2) are satisfied. 



4 Results, discussion, and 2D analogue 



As discussed in Section [T| we are interested in computing an upper bound (1.3) for the dimension 



Ajnin of the leading scalar in the OPE 0^ x 0^, universal for all unitary 4D CFTs. In [7], we have 
computed such a bound in the interval 1 < d < 1.3, using the sum rule of Section [2l2] truncated 
to the A^ = 6 derivative order. That bound is reproduced in Fig. [3j 

We now present the results of our latest study, obtained for larger values of A^. These result^ 
are plotted in Fig |6] as a collection of curves fwid), N = 6 ... 18, where the index A^ denotes the 
number of derivatives used to obtain the bound. The bound naturally gets stronger as A^ increases 
(see below), and thus the lowest curve fisid) is the strongest bound to date. In the considered 
interval 1 < d < 1.7 this bound is well approximated (within 0.5%) by 

/i8(rf) -2 + 0.77^/^ + 2. 17 + 0.437^/^ 7 = d-l. (4.1) 



To obtain the bounds of Fig. [6| we used the positivity argument from [7], as reviewed in Section 
^ Namely, for points lying on the curves Amin = fN{d) we are able to find a linear functional of 
the form (3.1) satisfying the positivity property (3.2 )|^ The numerical procedure that we use to 



find these "positive functionals" is described in some detail in Appendix A. 



^See Appendix B for the same results in tabular form. 

''Thus actually the bound is strict: Amin < fN{d), except at d ■ 
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Figure 6: Our main results. The solid curves are the hounds fN{d), N = 6... 18. 
The bounds get stronger as N increases, thus N = 6 is the weakest bound (highest 
curve), and N = 18 is the current best bound (lowest curve). The shaded region is 
thus excluded. The dashed curve foo{d) is an approximation to the best possible bound, 
obtained by extrapolating N ^ oo. The dotted line Amin = 2d is realized in a family 
of "generalized free scalar" CFTs, and is compatible with our bounds. 
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Figure 7: Same as Fig. |6| but with anomalous dimensions d — 1, A — 2 in logarithmic 
scale. The shaded region is excluded. 
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Several comments are in order here. 



We have actually computed the bound only for a discrete number of d values, shown as points 
in Fig. [6] The tables of these computed values are given in Appendix B. Behavior for d ^ 1 
can be better appreciated from the logarithmic-scale plot in Fig. [7j 

We do not see any significant indication which could suggest that the curves fN^d) do 
not interpolate smoothly in between the computed points. Small irregularities in the slope 
are however visible at several points in Figs. [GfT} These irregularities are understood; they 



originate from the necessity to discretize the infinite system of inequalities (3.2), see Appendix 
A for a discussion. In our computations the discretization step was chosen so that these 
irregularities are typically much smaller than the improvement of the bound that one gets 
ioT N ^ N + 2. 

For each the bound fN{d) is near-optimal, in the sense that no positive functional involving 
derivatives up to order A^ exists for 

A,,i,-2<(l-£)[/^(rf)-2]. 

We estimate e ~ 1% from the analysis of residuals in the fit of fN{d) by a smooth curve like 



in (4.1). 



On the other hand, by increasing A^ we are allowing more general functionals, and thus 
the bound fN{d) can and does get stronger. This is intuitively clear since for larger A^ the 
Taylor-expanded sum rule includes more and more constraints. 

Compared to the results of [7], the bound on the anomalous dimension Amin — 2 is improved 
by ~ 30 -i- 50% in the range 1 < ci < 1.7 that we explored. 

3. We have pushed our analysis to such large values of A^ in the hope of seeing that the bound 
saturates as A^ — oo. Indeed, we do observe signs of convergence in Figs. [6|7[ especially at 
d> 1.1. In fact, we have observed that the bounds fN^d) starting from N = 8 follow rather 
closely the asymptotic behavior 

fN{d)^Ud) + ^, (l<rf<1.7). 

An approximation to the optimal bound foo{d) can thus be found by performing for each d 
a fit to this formula. This approximation is shown by a dashed line in Fig. [6j From this 
rough analysis we conclude that the optimal bound on the anomalous dimension Amin — 2 is 
probably within ~ 10% from our current bound. 

4. We have fN^d) — > 2 continuously as d ^ 1. The point d = 1, Amm = 2 corresponds to the 
free scalar theory. 

We don't know of any unitary CFTs that saturate our bound at > 1, see the discussion 
in Section 6 of [7j. We know however a family of unitary 4D CFTs in which A^[^ = 2d and 
which are consistent with our bound (the red dotted line in Fig. [6]). This "generalized free 
scalar" theory is defined for a fixed d by specifying the 2-point function 



(0(x)0(O)) = \x 



-2d 



10 



and defining all other correlators of via Wick's theorem. This simple procedure gives a 
well-defined CFT, unitary as long as d > 1, which can be described by a nonlocal action 

Soc J d^x(t){d'')U. 

The full operator content of this theory can be recovered by studying the OPE (j) x (p. In 
particular, the leading scalar in this OPE has dimension 2d^ 



4.1 2D analogue 

Although our main interest is in the 4D CFTs, our methods allow a parallel treatment of the 2D 
case. The main characteristics of the 2D situation were described in Section 6.1 of [7j, here will 
briefly review them. 

1. At present we can only take advantage of the finite-dimensional SL{2, C) symmetry and not 
of the full Virasoro algebra of the 2D CFTs. In particular, our results are independent of the 
2D central charge c. 

2. The unitarity bounds for SL{2, C) primarie^in 2D have the form 

A>/, / = 0,1,2..., 

where / is the Lorentz spin. 

3. The SL{2,C) conformal blocks in 2D are known explicitly jTU] : 

9A,ii^^^) = ^ [/a+/(^)/a-K^) + ^)] • (4.2) 



Using the unitarity bounds, the known conformal blocks, and the sum rule (2.8), valid in 
any dimension, we can try to answer the same question as in 4D. Namely, for a SL{2, C) scalar 
primary of dimension d, what is an upper bound on the dimension Amin of the first scalar 



operator appearing in the OPE 0x0? I.e. we want a 2D analogue of the bound (1.3). Since the 
free scalar is dimensionless in 2D, the region of interest is (i > 0. 

Fig. [s] summarizes our current knowledge of this boundj^ 



The dotted line is the old N = 2 bound presented in [7]. The solid line is the iV = 12 
improved bound obtained by us{^ An numerical fit to this bound is given by: 

A.Sd + Sd"^ -87d^ + 2300d\ rf< 0.122, 
0.64 + 2.87ci, d> 0.122. 



^This theory can also be realized holographically by considering a free scalar field of a particular d-dependent 
mass in the AdS geometry and taking the limit in which 5D gravity is decoupled. We are grateful to Kyriakos 
Papadodimas for discussions about the generalized free scalar CFT. 

^Known as quasi-primaries in 2D CFT literature. 
^^See Appendix B for the results in tabular form. 

^"'^We are grateful to Erik Tonni for providing us with the large A, / asymptotics of the 2D conformal block 
expansion coefficients, necessary to obtain this bound. 
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Clearly, the improvement compared to [7] is significant. 

It is interesting to note that in 2D we have observed a much faster convergence for increas- 
ing than in 4D. In fact, aheady with = 6 it is possible to obtain a bound rather close 
to the one shown in Fig. [8} although with a shghtly rounded "knee" . We have also computed 
several points for = 16 and haven't seen much improvement. 

• The dashed line and scattered crosses correspond to various OPEs realized in explicit exam- 
ples of exactly solvable unitary 2D CFTs (minimal models and the free scalar theory), see 
[7j. They all respect our bound. 



2.0 
1.5 
J 1.0 
0.5 
0.0 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 

d 

Figure 8: See the text for an explanation. The red cross denotes the position of the 
Ising model, the black crosses marked ijj, ip"^ correspond to the OPEs realized in the 
higher minimal models, as in Fig. 15 of The shaded region is excluded. 




It is instructive to compare this plot with its 4D counterpart. Fig. [6j While we do not know of 
any CFTs saturating the 4D bound, the 2D unitary minimal models A1(m, m -|- 1), m = 3, 4, . . ., 
contain the OPEs 

V^x V = l + V^' + ..., A^ = ^- ^ A^2=2 (4.3) 

2 2(m + 1) m + \ 

which come quite close to saturating the 2D bound. 

More precisely, our 2D bound starts at (0, 0) tangentially to the line A = Ad realized in the 
free scalar theory, then grows monotonically and passes remarkably closely above the Ising model 
point (Act, Ag) = (1/8, 1). After a "knee" at the Ising point, the bound continues to grow linearly. 



passing in the vicinity of the higher minimal model points (4.3) 
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It is curious to note that if we did not know beforehand about the Ising model, we could have 
conjectured its field dimensions and the basic OPE a x a = 1 + e based on the singular behavior 
of the 2D bound at d = 1/8. 

On the other hand, nothing special happens with the 2D bound at the higher minimal model 



points, it just interpolates linearly in between Most likely, this does not mean that there exist 



other unitary CFTs with intermediate operator dimensions. Rather, this behavior suggests that 
the single conformal bootstrap equation used to derive the bound is not powerful enough to fully 
constrain a CFT. 

In comparison, it is a bit unfortunate that the 4D bound does not exhibit any singular points 
which would immediately stand out as CFT candidates. Nevertheless, if we assume that the shape 
of the 4D bound is a result of an interpolation between existing CFTs (as it is the case in 2D), 
we may conjecture that the upward convex behavior of the functions fNid) in Fig. [6] is due to the 
presence of a family of points satisfying the sum rule that can correspond to exact CFTs. This 
observation, though speculative, shows how the presented method can provide a guideline in the 
study of 4D CFTs. 



5 Future research directions 

The results of this paper and of [7] open up many interesting research directions, which we would 
like to list here. 

First, there are several important problems in 4D Conformal Field Theory which can be ana- 
lyzed by our method and its simple modifications. For example: 

1. One should be able to derive a generalization of our bounds in the situation when the CFT 
has a global symmetry, and we are interested in the lowest dimension singlet appearing in 
the OPE. This is going to have phenomenological implications by constraining the so-called 
conformal technicolor scenarios of Electro Weak Symmetry Breaking [llj. This connection 
was extensively discussed in [7]. 

2. One should be able to derive model- independent bounds on the size of OPE coefficients. This 
is going to be relevant for discussions of 'unparticle self-interactions' [12j, in the context of 
unparticle physics scenarios 



Second, the method can also be used in 2D Conformal Field Theory, as was already demon- 



strated in Section 4A The main interest here lies in potential applications to string theory. We 
will now briefly describe two such applications. 

Physical states of (super) string theory are in 1-1 correspondence with Virasoro primary oper- 
ators of a 2D CFT living on the string worldsheet. The mass of a string state (in string units) is 



^^The straight line fitting the bound would cross the dashed Free theory line just above d — 0.5, which is the 
accumulation point of the minimal models A4{'m,m+ 1). For larger values of d we expect that the bound modifies 
its slope and eventually asymptotes to the Free line. 
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related to the corresponding primary operator dimension A via 

= A - 2 . 



We are considering closed string theory for concreteness. When strings propagate in fiat space, 
the CFT is solvable and the full spectrum of operator dimensions is known. Realistic string 
constructions require compactifications of extra dimensions. In some examples, such as toroidal 
compactifications, the CFT is still solvable. In others, such as superstring compactifications on 
a generic Calabi-Yau three-fold, the CFT cannot be solved exactly. All what is generally known 
is the spectrum of the massless states, which can be obtained in the supergravity approximation. 
Of course we expect the massive string states to be always present, but just how heavy can 
they be? We know from the experience with toroidal compactifications that it is impossible to 
completely decouple the massive states: as the compactification radius R ^ 0, the Kaluza-Klein 
states become more massive, but the winding modes come down. Clearly, massive string states 
are crucial for the consistency of the theory. What exactly are they doing? A partial answer may 
be that without their presence, 4-point functions of the massless state vertex operators would not 
be crossing-symmetric. If this intuition is right, it could be used to obtain model-independent 
bounds on the lightest massive states in string compactifications, generalizing the well-known 
bounds valid for toroidal compactifications. A similar in spirit general prediction of string gravity, 
although in a different context and by using different methods, was obtained recently in |14j . 

When working towards results of this kind, it may be necessary to generalize our methods so 
that the information about the 2D CFT central charge, which is fixed in string theory, can be 
taken into account. In practice, one needs an efficient method to evaluate the full Virasoro con- 
formal blocks. While no closed-form expression as simple as Eq. (4.2) is known, Zamolodchikov's 
expansion (see [6j) can probably be applied. 

Finally, as mentioned above in the 4D context, it should be possible to derive model-independent 
bounds on the OPE coefficients. Such results must be accessible via a simple modification of our 
method, in particular the full Virasoro conformal blocks are not needed here. One can then apply 
such bounds to the dimension 2 operators corresponding to the massless string states (in an arbi- 
trary compactification). Via the usual dictionary, this would then translate into general bounds 
on the tree-level coupling constants in the low-energy string effective actions. 



6 Summary 

Prime principles of Conformal Field Theory, such as unitarity, OPE, and conformal block decom- 
position, imply the existence of an upper bound f{d) on the dimension Amin of the leading scalar 
operator in the OPE (p x <f), which depends only on (p^s dimension d. 

Moreover, there is an efficient method which allows numerical determination of f{d) with 
arbitrary desired accuracy. The method is based on the sum rule, a function-space identity 
satisfied by the conformal block decomposition of the 4-point function {(j)(f)(f)(j)) , which follows from 
the crossing symmetry constraints. In practical application of the method the sum rule is Taylor- 
expanded: replaced by finitely many equations for the derivatives up to a certain order N. The 
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bound f{d) improves monotonically as more and more derivatives are included. In [7], where the 
above paradigm was first developed, we numerically computed the bound for = 6. 

The present paper extended the study of [7] to higher N. The goals were to improve the bound, 
and perhaps to approach the best-possible bound in case a convergence of the bound is observed. 

Our analysis went up to = 18, see Fig. [6| and we have achieved both goals. First, in the 
range 1 < d < 1.7 that we explored, the bound on the anomalous dimension Amm — 2 is improved 
by 30 -7- 50% compared to the results of [7j. Second, we do observe signs of convergence of the 
bound. We believe that our current results are close (within ~ 10%) to the best ones achievable 
with this method. 

The results of this paper and of [7] suggest several interesting research directions, connected 
with phenomenology and, via the 2D analogue of our method, with string theory (see Section Isl). 
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A Details about numerical algorithms 

We now discuss in more detail the issues introduced in Section |3| namely how we can find in 



practice a linear functional A[F] of the form (3.1) satisfying the positivity property (3.2). We will 



first describe the general procedure and how it can implemented in a computer code, and then 
mention possible algorithmic improvements and shortcuts that we found useful in our analysis. 

Given the complexity of the functions Fd,A,i, the search for a positive functional is too hard a 
task to be attacked analytically. As already mentioned, we reduce the complexity of the problem 
by looking for a functional which is a linear combination of derivatives up to a given order A^. 
The derivative are taken w.r.t. the selfdual point X = Y = 0, since the sum rule is expected to 
converge fastest around this point and, in addition, the functions F^ a.K^'^) even in both 



arguments. The choice of the functional (3.1) simplifies our task enormously since we can now 
work in a finite dimensional space, and the only information concerning Fd^A,i that we need are 
their derivatives up to a certain order. Put another way, the F-functions are now considered as 
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elements not of a function space but of a finite-dimensional vector space W, s = N{N + 6)/8. 



The sum rule (2.8) in this picture represents a constraints on these vectors that, in any CFT, 
must sum to zero. This interpretation is discussed in details in [7]. Here we adopt an equivalent 
point of view in terms of the dual space of linear functionals defined on M'^ since we find this 
prospective closer to the method used to obtain numerically A[F]. 

Let us fix the notation. We define the s-dimensional vector of Taylor coefficients: 

J^o[d,6,l] = l^—Fj,"^':^\m,neven,2<m + n<N\ , (A.l) 
[m!n! ' ' J 

^dA'f' = 9xdYFd,A,i\x=Y=oj 5 = A — / — 2, 
and the same vector normalized to the unit length: 

J^[d,6,l]^r^ , (A.2) 
Ik oil 

where the norm \\J-'o\\ is the usual Euclidean length of the vector jFg. 

We form the vector JFq out of the Taylor coefficients of the function i^d,A„z rather then of 
its derivatives, because this way all elements turn out to have approximately the same order of 
magnitude, which is preferable in the subsequent numerical computation. The definition of the 
normalized vector JF serves the same purpose. Indeed, as explained in the following, our numerical 
analysis consist in finding a solution of a system of linear inequalities where the coefficient are 
given by the elements of J-'[d,6,l]. The solution is more accurate and easier to extract if all the 
coefficient are of the same order of magnitude. Since the existence of the functional A is not 



affected by these rescalings, we opted for the definition (A.l), (A.2). 

According to the positivity property we look for a functional which is strictly positive on all 
but finitely many vectors J-'[d,6,l]. Let us fix the dimension d of the scalar (pd- Then each pair 
A , / identifies the semi-space of [W)* of the functionals positive-definite on the vectors J-'[d, 6, /]; 



let us call this open sets Ud^^i- With this notation the positivity property (3.2) can be restated 
in the following way: If for fixed d and A^in 



n 



Uj,A., 9 , (A.3) 



A>i+2,«=2,4,... 

then the sum rule cannot be satisfied. The issue is thus to be able to check whether the intersection 



(A.3) is non-empty, and to compute the smallest Aynm for which this is the case. 

Clearly it is not possible nor needed to check all the values of A as required by the condition 



(A.3). We can indeed consider only a finite number of them and check if they admit the existence 
of a functional or not. This can be achieved with a double simplification. First, we consider values 
of A and / only up to a given maximum value ("truncation"), and secondly, we discretize the kept 
range of A ("discretization"). The truncation does not produce a loss of information since we 
take into account the large A and / contributions using the asymptotic expressions computed in 
Appendix D of The discretization step requires special care, see below. 
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We used Mathematica 7 to perform the computations. The algorithm to extract the smallest 
value of Amin proceeds in several steps: 

1. Setting up an efficient procedure to compute vectors J^[d, 5, /]. 



2. Selection of the /'s and 5's to be used in checking the positivity property (|A.3|). For concrete- 
ness we report here the range of /, 5 that we were including: 
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2 < / < /max = 50 : < 5 < 200 , (A.4) 
/ = : = Amin - 2 < 5 < 200. 



For each / the range of 5 was discretized, and a discrete set of points was chosen, called F^ 
below. The derivatives of the F-functions approach zero as 5 — > cxd and reach the asymptotic 
behavior for sufficiently large values. We take a finer discretization where the function are 
significantly varying while we can allow to increase the step in the asymptotic region. More 
details are given below. 

Reduction to a Linear Programming problem. With only a finite number of equations to 
check, the determination of the intersection of the Ud,/\,i becomes a standard problem of 
Linear Programming which can be solved in finite amount of time. Hence we look for a 
solution of the linear system of inequalities 

K[T[d, 6, 1]] = J2 Am.n^m,™ > , (A.5) 
5 eTi, 1 = 0... /„jax • 



Clearly, the coefficients Xm,n are related to those appearing in (3.1) by a trivial rescaling 
depending on m,n: 

Am,n 1^^'1^^'Xm,n 

Further, the asymptotic behavior of the F-functions (see below) tells us that for large 6 the 
inequality is dominated by the (A^, 0) derivative 

A[FdA,i] ^ A^.oi^iS? {S » / » 1), (A.6) 
hence Aat^q needs to be positive. By an overall rescaling of A we can always achieve 

Aiv,o = l, (A.7) 
which we choose ClS db normalization condition. 

Extraction of the smallest Amin for which a positive functional exists. We begin by selecting 



two points 5min = ^1 and 5min = ^2 > ^1 (see (A.4)) such that we know a priori that in the first 



case a positive functional does not exists, while in the second case it does Starting from 



^^In some cases 2 < I < 200 was needed to obtain a functional which would later pass the positivity check on 
the non-included values of I, see below. 

^''We can choose these points blindly as Si = 0, S2 ^ 1, however prior experience can suggest a choice closer to 
the final Smin 
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these values we apply the bisection method to determine the critical 5yam up to the desired 
precision: we test if a functional exists for 6mm = (^2 + 5i)/2 and we increase or decrease the 
extremes of the interval [5i,52] depending on the outcome. The procedure we follow is such 
that in the end the critical Smm is contained in an interval of relative width 10~^, i.e. we 
terminate if ^2 — ^ 10~^5i. The plots and the tables presented in this work correspond to 
the upper end 62 of the final interval, i.e. to the the end for which we have found a functional. 

Let us now come back to the point 1. Although computation of the derivatives can be carried on 
by brute force Taylor-expanding the F-functions, we can save time decomposing the computation 



in various blocks. From equation (2.8) we see the rather simple dependence on the parameter d, 



which translates in a polynomial dependence once the function is Taylor-expanded in X and Y. 
We therefore separately computed the dependence on d once and for all as a matrix M{d)mn\ij- To 
compute Taylor coefficients of F-functions, this matrix is contracted with two vectors containing 
one-dimensional Taylor coefficients of the function kp{x), see (2.4). The latter derivatives are 



pre-computed for several values of (3 with a fine step and stored. For definiteness we report the 
interval we used: 

< /5 < 10^ step: 10-3 . (A.8) 

For larger j3 we made use of the analytic expression of the asymptotics instead of computing the 
derivatives numerically (see below). 

Finally let us discuss the choice of the discretization and the truncation in A and /. This step 
is of fundamental importance in order to reduce the time needed to perform the computation. 

In Appendix D of \J\ it is shown that for large values of 5 and I the functions -F^^a,; approach 



an asymptotic behavior. We have checked that outside the range of values (A. 4) we can safely use 
the approximate expression 



-,(m,n) 



const. {2\f2) 



m+n+2 



(m + l)(n + 1) ' 

For large /, 6 the vector J^[d, S, I] is dominated by the components where m + n assumes the highest 
allowed value N. Hence we can take into account this large /, 6 behavior imposing additional 
constraints: 



A[J-e] > , 



(cos e+sin 6')""+l cos 8"+^ 
(m+l){n+l) 





if 



m + n = N 
otherwise 



tan 6 



6 

r 



where we have dropped irrelevant positive constants not depending on m , 



n. 



Now comes the discretization: in the range of values (A. 4), as well as in the interval 6 G [0, 7r/2], 
we can allow to take only a discrete, finite number of points. For 6 we take a fixed small step. 
However, for 6 we try to concentrate the points in the region where the unit vector J-'[d, 6, 1] is 
significantly varying. A measure of this is given by the norm of its derivative w.r.t. 6: 
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AT 



_d_ 
06 



J^[d, S, I] 
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In practice the derivative d/dd is evaluated by using the finite-difference approximation. 
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We discretize by taking the spacing between two consecutive values of S equal c/Af, where c is a 
small fixed number (c = 0.02 -=-0.05 was typically taken in our work). Clearly when the unit vector 
is slowly varying the discretization step is large, while it is refined where it is changing rapidly, 
and where presumably more information is encoded. Typically we get about a hundred 6 values 
for each I, but only a few dozen of those above 6 > 50. 

The sets Ti, one for each /, of values of S obtained in this way are the ones referred to at 
point 2 above. In constructing the linear system that we use at point 3 we consider additional 
intermediate points between two subsequent (5's. In order to understand why we do this, let us 
assume that we have found a functional A which is positive for all the values of 6 contained in Ti. 
Since we considered a discrete set of values, it may and actually does happen that for intermediate 
values of 6 (which were not included in F/) the functional becomes slightly negative. In [7] this 
issue was solved looking for solution of the form A[F(;i a,/] > so that for intermediate values this 
condition could be violated but the positivity was safe. In the current work we found it more 
convenient to build the linear system in the following way: 

• for each 6 eTi = {6i, 5i, ....} we evaluate the vector J^[d,6,l]. 

• for any two consecutive points we consider the first-order Taylor expansion of the 
vector jF[(i, 5, /] around 6 = 6i and evaluate it at half-spacing between 6i and 

J'y2[d,S„l]=J'[d,S„l] + ^J'[d,5,,l] (A.9) 



and we add the constraints A[jFi/2] > to the linear system (A. 5). 



These additional constraints become important to keep the functional positive near the (5's for 
which the inequalities A[jF] > are close to saturation, while they are redundant away from those 
points. Indeed, assume that for some 6i and the functional is exactly vanishing. Then at the 
intermediate point the functional would be strictly negative, which is not allowed. However, in 
the presence of the additional constraint A[jFi/2] > this cannot happen, since A[JF] is generically 
a convex function of 6 near the minimum. See Figure |9] for an illustration. Thus we can be certain 
that the found functional will be positive also for those 6 which were not included into F^. 

This certainty has a price. Namely, the opposite side of the coin is that the added J-'i/2 
constraints are somewhat stronger than needed, and the bigger the discretization parameter c, 
the bigger the difference. As a result, the found critical value of Amin will be somewhat above 
the optimal critical value, corresponding to c — 0. This observation explains why the curves 
in Figs. [6|7l have small irregularities in the slope. These irregularities could be decreased by 



decreasing the value of c. 

Several comments concerning the numerical accuracy are in order. The components of the 
vector J^[d,S,l] have been computed using standard double-precision arithmetic (16 digits). As a 
consequence all the numerical results must be rounded to this precision. In particular, quantities 
smaller that 10~^^ are considered zero. 

In addition, the built-in Mathematica 7 function LinearPrograjiiming, which we used, has 
an undocumented Tolerance parameter. Most of the computations were done with Tolerance 
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A[F[d,6,l]] 



A[F[d,6,l]] 




Figure 9: Imposing the positivity of the functional on a discrete set of points, it could 
happens that the intermediate points don't satisfy A[J-'[d, 6, 1]] > (on the left). However 
adding the constraint A[J^i/2[d,6i,l]] > 0, see (A. 9), we can be sure that the functional 
is positive on all the neglected points (on the right). 



equal 10~^ (default value). However for = 16 and = 18, and for d < 1.1, we found that 
LinearProgramming terminates prematurely, concluding that no positive linear functional exists, 
even for some values of ^min for which a positive functional for smaller was in fact found. The 
problem disappeared once we set Tolerance to a lower value (10~^^). In our opinion. Tolerance is 
probably the so-called pivot tolerance, the minimal absolute value of a number in the pivot column 
of the Simplex Method to be considered nonzero. Recall that a nonzero (actually negative) pivot 
element is necessary in each step of the Simplex Method [12] . This interpretation explains why 
the above problem could occur, and why it could be overcome by lowering Tolerance. 

As described above, our numerical procedure has been designed to be robust with respect to 
the effects of truncation and discretization. In addition, for each d, we have tested the last found 
functional (i.e. for (5mm at the upper end 62 of the final interval [61, 62]) on the much bigger set of 
6,1: 

2</< 500, < 5 < 500 ,step = 0.1, 

^ = , 5min < 5 < 500 , step = 0.1 , (A.IO) 

and found that indeed A[jF] > 0, within the declared 10~^^ accuracy. 

Finally, we have checked that in all cases the found functionals A are such that the inequality 
A[jF] > is in fact strict: A[jF] > 0, for all but finitely many values of S and /. Thus they satisfy 
the requirements stated in Section [3] 



B Tables 



Table 1 contains the sequence of 4D bounds fN{d),N = 6...18, for a discrete set of points in the 
interval 1 < d < 1.7. Table 2 contains the 2D bound fi2^\d) for a discrete set of points in the 
interval < d < 0.35. Figs. |6| [7], |8]are based on these tables. 
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A text file with the unrounded versions of Tables 1,2 and the functionals used to derive these 
bounds is included in the source file of this arXiv submission. The Mathematica codes can be 
obtained from the authors upon request. 



fN{d) - 2 



d 


- 1 


= 6 


N ^8 


= 10 


A^ = 12 


A" = 14 


A^ = 16 


A^ = 18 





01 


0.2045 


0.2025 


0.1385 


0.1356 


0.1006 


0.08955 


0. 





02 


0.3055 


0.2956 


0.2086 


0.1965 


0.1636 


0.1485 


0.1425 





03 


0.3895 


0.3676 


0.2636 


0.2496 


0.2126 


0.1956 


0.1916 





04 


0.4646 


0.4315 


0.3165 


0.3006 


0.2506 


0.2366 


0.2286 





05 


0.5358 


0.4878 


0.3665 


0.3467 


0.2896 


0.2786 


0.2596 





07 


0.6688 


0.598 


0.4576 


0.4296 


0.3927 


0.3627 


0.3397 





1 


0.848 


0.7348 


0.5847 


0.5486 


0.5026 


0.4705 


0.4448 





15 


1.106 


0.957 


0.7789 


0. 


0.6807 


0. 


0. 





2 


1.354 


1.163 


0.9635 


0.9118 


0.8396 


0.8098 


0.7757 





25 


1.597 


1.362 


1.142 


0. 


1.002 


0. 


0. 





3 


1.841 


1.56 


1.319 


1.249 


1.163 


1.127 


1.085 





35 


2.085 


1.759 


1.494 


0. 


1.317 


0. 


0. 





4 


2.329 


1.957 


1.669 


1.578 


1.473 


1.436 


1.389 





45 


2.561 


2.159 


1.843 


0. 


1.63 


0. 


0. 





5 


2.785 


2.365 


2.046 


1.907 


1.786 


1.745 


1.693 





55 


2.994 


2.571 


2.198 


0. 


1.942 


0. 


0. 





6 


3.213 


2.78 


2.381 


2.238 


2.1 


2.054 


1.995 





65 


3.441 


2.995 


2.561 


0. 


2.258 


0. 


0. 





7 


3.638 


3.21 


2.744 


2.574 


0. 


2.366 


2.297 



Table 1: 4D results. The table contains anomalous dimensions, i.e. d—1 and fN{d) — 2 
are given. A zero entry means that the bound for this d and N has not been computed. 
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rl 


f 
J 


12 l"J 


0.01 





04254 


0.02 





08752 


0.03 





1356 


0.04 





1865 


0.05 





2406 


0.06 





2998 


0.07 





3676 


0.08 





4455 


0.075 





4055 


0.08 





4456 


0.085 





4907 


0.09 





5408 


0.095 





5947 


0.1 





6577 


O.iOo 





1 6 


0.11 





8037 


0.115 





8716 


0.12 





9418 


0.123 





9803 


0.125 


1 


001 


0.127 


1 


008 


0.13 


1 


017 


0.15 


1 


073 


0.2 


1 


214 


0.25 


1 


357 


0.35 


1 


647 



Table 2: 2D results. 
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